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ABSTRACT 

A coupled system of two generalized logistic maps is studied. In par- 
ticular influence of the coupUng to the behaviour of the Julia set in 
two dimensional complex space is analyzed. It is proved analytically 
that the Julia set disappears from the complex plane uniformly as a pa- 
rameter, which interpolates between chaotic phase and integrable phase, 
approaches to the critical value as far as the coupling strength satisfies 
a certain condition. 
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I. INTRODUCTION 



The study of nonlinear dynamical systems in the last two decades uncovered the existence of hierarchy of 
completely integrable models in one hand and of universal structure of nonintegrable models on the other 
hand. The former is characterized by concepts such as solitons, inverse methods, Yang-Baxter equations, 
infinite dimensional symmetries etc., whereas the terms such as chaos, fractals, strange attractors, etc. 
characterize the latter. 

Since there exists no general method to analyze nonlinear equations the completely integrable models 
are quite exceptional among other nonlinear systems. Nevertheless they appear commonly in fundamental 
theories of physics, such as Einstein's theory of gravity, string theory of elementary particles, solvable 
lattice models in statistical physics and soliton models in various fields. They are easily distinguished 
from the behaviour of chaos solutions of nonintegrable systems. They must be also distinguished from 
stable solutions of nonintegrable systems. 

A Julia set is a set of points of initial values on the complex plane of dependent dynamical variables 
whose iterational mapping never converge^ A chaos will be observed experimentally if a Julia set 
intersects the real axis in the complex plane. Physical phenomena are expected being stable if the Julia 
set stays far from the real axis. The existence of a Julia set somewhere on the complex plane prevents 
the system to be solved analytically, besides some exceptional cases. Conversely there is no room for a 
Julia set to appear in the complex plane if the system is completely deterministic and predictable. From 
this view point the Julia set is an invariant of the dynamics, hence is a proper object which characterizes 
complete integrable systems among others. 

Some topological aspects of Julia set have been studied for simple cases from mathematical view points 
1^. It must possess full of information of the dynamics. Their physical significance is, however, not 
known. We like to emphasize here that the Julia set is very sensitive topological object against to any 
small difference of parameters. For instance two systems which are in chaotic orbits independently, often 
turn to a stable motion when they are coupled weakly. This can be understood that the Julia set is 
removed from the real axis as the coupling increases from zero. 

In the previous papers , |Q we studied the behaviour of the Julia set at the point of transition from 
a nonintegrable phase to an integrable phase by using a simple model. It includes both the integrable 
and nonintegrable logistic maps as special cases of a parameter which interpolates them. This parameter 



^This definition of a Julia set applies in one dimensional maps. Since there has been not known unique definition 
of Julia set in higher dimensional maps, what we call here a Julia set is somehow vague. More precise definition 
appropriate to our purpose will be given in §4. 
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enables us to study the system analytically near the critical point of transition. It was shown that the 
Julia set which was a Cantor set concentrated along the real axis became a countable number of discrete 
points as a parameter approached to the critical point of the transition Q. Moreover in |^ we have 
shown that the Julia set of the generalized logistic map (GLM) converges uniformly to a certain regular 
orbit in the integrable limit. 

The purpose of this paper is to study the behaviour of the Julia set under the influence of the coupling. 
To be specific we consider the generalized logistic map (GLM) as a unit of dynamical system. Our 
strategy of studying the coupling dependence of the Julia set is to find a model of coupling as simple as 
possible, so that we can see analytically the influence of the coupling. Therefore we study in this paper 
a system in which there are two coupled GLM's. 

When two units of GLM's are coupled, the Julia sets appear on each complex plane of the variables. 
Since they are correlated with each other through the coupling we are to study their behaviour in the 
hyper plane of two complex variables, or four real variables. 

In spite of the complexity arose from the hyper complex plane, we will show that our coupling model 
of two GLM's enables us to extract certain informations about the Julia set. In particular we will show 
that, in addition to those in the case of single map, there appear new fixed points of the map which cause 
a deformation of the Julia set. We will show that we can study such deformation analytically to some 
extents. We also prove that the Julia set converges uniformly to the orbits of integrable map, as far as 
the value of the coupling strength satisfies a certain constraint. 

Our paper is organized as follows. We will review some results of the case of single GLM in the 
next section. In §3 we introduce the coupled system and study its basic features. The main results are 
discussed in §4, where the coupling dependence of the Julia set is analysed in detail. In order to supply 
informations we will present numerical results of some projections of the Julia set, so that we can see 
visually the coupling dependence. The last section is devoted to discussions. 
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II. JULIA SET OF SINGLE GENERALIZED LOGISTIC MAP 



The Julia set itself does not depend on initial values, but is uniquely determined by dynamical equation. 
Since it is a highly complicated set of points in the complex plane analytical study of its behaviour has 
been limited to simple models Q. In particular little is known about properties of Julia set in the case 
of coupled systems. Therefore, before studying the coupled case, it will be worthwhile to summarize 
properties of the Julia set of a single generalized logistic map (GLM), which we investigated in our 
previous work 

A. General feature of the Julia set of single GLM 

In the previous paper jl| , we studied a map which reduces to the logistic equation in the continuum 
limit and interpolates integrable and nonintegrable logistic maps: 

zi+i = fizi):=ii^-^-^j^j—^^; zieC, leZ. (1) 

Here /x, i/ and 7 are complex parameters. When 7 = 1 this is the standard logistic map. When 7 = it 
becomes a Mobius map, hence we called it the integrable logistic map (ILM). The latter version of the 
logistic map was first studied by Morisita as an alternative difference analogue of the logistic equation. 
It has a unique solution given by 

zi = — (2) 

1 + V- zo 

1 

for an arbitrary initial value zq G C . The continuous time limit of (|^) can be taken by replacing 

zi = — ■ u{t), hm = — 37"W' 

v + ^fi — h^a n u + 7/i — 71/ at 

jj, — I + ah (3) 

in the limit of /i — s- 0, from which follows the logistic equation : 

— =au{l-u). (4) 

The merit of studying the map is that we can investigate analytically the transition of the dynamical 
system from nonintegrable phase to the integrable phase. 

We will use more convenient form, in the place of (Q), the standard parameterization of the rational 
map of rank 2 

n^) = r^o/o0(z) = i(i±^e^« (5) 



A = ^e-^ A' = n_J^^ln±I^e^o^ (6) 

after the proper transformation : 

^^^^ ^ P)x 

(1^7 — V — nj)x + (i^7 — v — 7)/ze^''^ 

In this form of the map ILM becomes extremely simple : 

F{z) = fxz, if 7 = 0. (8) 

In the new parameters the integrable limit 7 = corresponds to A A' = 1. 

The most fundamental analytical objects are fixed points. The GLM has three fixed points, 0, 00 and 

Zp 

_ A - e-'^ 

whose properties can be easily deduced from the values of derivative at those points. When its absolute 
value is less than 1, the point is attractive and when it is greater than 1, the point is repulsive. 

As it was explained in the introduction, the Julia set is an invariant of the mapping. Therefore it 
should provide us useful knowledge about global nature of the system. 

A Julia set is a set of points of initial values whose iterational maps neither converge nor diverge. There 
are several equivalent methods to calculate the points. We will adopt the following definition of the Julia 
setQ: 

Definition : If ip is a map of z into (p{z), and exists its inverse map ip~^ , then the Julia set J{ip) 
of if is defined by the following set of points : 

J{^) |z| y ^-\zj),^zj e J(^)| . (10) 

If zj is a point of the Julia set, its every image belongs to the Julia set. Conversely, if points in a 
neighborhood of the Julia set do not belong to the Julia set, then all images spread to whole complex 
plane. Therefore the definition makes sense. 

From our definition of Julia set described above, we can create the Julia set by considering its inverse 
mapping. The inverse mapping of GLM is given by 



zi = F~\zi+i) 

^ ^{pzi+, - A) ± ^\/{pzi+, - A)2 + 4zi+ie-»», (11) 



**This is not the conventional definition of the Julia set. Nevertheless we adopt this definition which is 
natural and more convenient to understand the transition to the integrable pheise. 



5 



where we defined 

p:=A'e-'^. (12) 

Attractors of this mapping constitute the JuUa set. The feature of the mapping is determined by dynamics 
of the system. 

It is well known that the behavior of a dynamical system, hence the nature of the Julia set, is sensitive 
to the values of the critical points, the zeros of the derivative of the mapping. In our problem, by solving 
F'{z) = 0, they are 

zcrt = ^{-l±Vl-XX'). (13) 

The mapping is discriminated by the value of 

e:=A'A-l. (14) 

To be specific we consider the case of real parameters, i.e.. A, A' are real and set 9 = 0. Under these 
conditions, e is real, and there arise three types of critical points, depending the values of e: 

1. Two complex critical points (e > 0) 

There is no critical point on the real axis and graph becomes monotonically increasing two curves. 
The inverse mapping transforms a real value to a real one. If A' > 0, one of the fixed points 0, which 
is on the real axis, is a repulsive one. Hence the whole iterational images of this point belong to 
the Julia set which are constrained on the real axis. In fact it becomes a Cantor set as mentioned 
in ref. . 

2. One critical point (e = 0) 

Two critical points join into one at the singular point z = —A, and dynamics becomes linear. In 
this case the Julia set has gone and system becomes integrable. We will call e the integrability 
parameter in what follows. 

3. Two real critical points (e < 0) 

The graph has two peaks along the real axis and they can be considered as two independent 
monotonic systems. In this case one of the critical points in between two fixed points and Zp is 
important. Many patterns of the Julia set appear in the complex plane and exhibit complicated 
structures such as Jordon curves, Cantor sets or dentrites, depending on location of the critical 
points. 
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B. Integrable limit of the Julia set of single GLM 



In our previous paper |g] we have shown in the case of single GLM that the Juha set converges uniformly 
to the set of points of the Mobius maps in the integrable limit. We will summarize briefly the proof in 
this subsection. 

Let us first consider the integrable case, i.e., e = 0. If |A| > 1 the origin is a repeller. The inverse map 
( |ll|) started from z = turns to be 



zi = 



which yields an orbit consisting of the set of points 

J^^*^ = {-p"A n e N}. (16) 

We now consider the cases in which e does not vanish. The inverse map ( |ll| ) becomes multivalued and 
much complicated. Nevertheless we could prove the following fact: 

All points of the Julia set approach uniformly to J^^*^ in the integrable limit. 

In order to see the behavior of the Julia set at small values of e, it is convenient to rewrite (111) as 



P''i^)^[ifJ) A: z^pz + E{z), B:z^-\^E{z) (17) 



where 



Ei.):=\y.^X)Ul-^i^-l]. (18) 



E{z) vanishes for small values of e. In fact we can show 

\E{z)\ < R,, (19) 

for all values of e. Here is given by|^ 

Re = ^^^VW\(VW\ + VV+M)- (20) 



lA'l 

which tends to zero like -^/jej as e approaches to zero. It is remarkable that i?e does not depend on the 
argument of the map, but is determined only by the parameters characterizing the map. 
We can prove our claim as follows: 



^This result, which we owe to Hisao Konuma, improves the one of our previous paper Q. 
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1. Using A and B defined in ( p^ ) tlie nth iteration of F ^ yields 



F-"(z) = {A^'B-'^A''' ■ --B^-iz)] vi + V2 + --- + Vn = n} ■ 



(21) 



2. For any X an element of the form A^BX, s = 0, 1, 2, • • • , n — 1 in ( |21| ) lies in the neighbourhood of 



\A'{BX)+p'\\ 



< 



-p'E{X) + p^E [A'-^-^BX] 
1-|P| 



fc=0 



s==0,l,2,---,n-l. 



(22) 



Since we assumed |A| > 1, \p\ = |A'| < 1 is satisfied as long as e is small. Hence the right hand side 
of (p^ ) is finite for all s. 

The exception is the case y4."(0), for which |y4"(0)| = holds since z = is a fixed point of the map 
A. 



3. From 1. and 2. we conclude 



1 



F-"(zo)-(-p"A) <- ^Re^O, ^n e^O 

In our present problem the Julia set is the collection of all F~"(0), i.e., 

oo 

jGLM ^ Q g jGLM^ 

n=0 



(23) 



(24) 



hence the claim is justified. 
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III. INTERFERENCE OF TWO GLM'S 



We now study the coupling of GLM. The most general linear coupling is given in the following form 

4\=E5..^^(4'^) (25) 

k 

with arbitrary constants {gjk} and F{z) being a mapping function. We should recall here the fact that 
the Julia set has been studied mainly from mathematical point of view, and restricted to the case of 
small number of variables. Since we are interested in studying influence of the coupling to the Julia set 
generated by the function F ^ we will restrict our discussion to the simplest case of the coupling, namely 
to the coupling of two units. Hence the system to be studied in the following is then 

zi+i = Lp{zi, z'l) = (1 - g)F{zi) + gF{zi), 

- ^{z'l, zi) = gF(zi) + (1 - g)F{z[) (26) 

A. Fundamental feature of the mapping 

First let us notice some results which can be drawn from this expression by a glance. It has an 
invariance under the change of the coupling strength g 1 — g if the variables zi and z[ are exchanged 
simultaneously. Therefore it is sufRcient to investigate the limitted region g < ^. When g — two 
equations of ( ^6|) coincide, hence the map reduces to the one of single GLM. 

( p6|) is a map from two dimensional complex space on itself in general. There exists, however, a special 
subset of maps which close among themselves and remain essentially one dimensional. It occurs once z; 
and z[ coincide each other. Then, as it is clear from the expression of ( |2^ , all subsequent mapping give 
the common values to them, and the mapping becomes equivalent to a single map. This phenomenon 
is not special only for (|2^), but it happens to all linear mapping of the form ( p5| ) under the constraint 
Y2=i9jk = 1, j = 1,2,- • - ,71. 

The most basic informations are given by fixed points of the mapping. In the case of coupled maps it 
is not clear at all if there exists a way to find them analytically. Therefore it is remarkable that in the 
case of two unit coupling of GLM we can find them analytically, as we see below. This fact enables us to 
study dynamics of the system in detail. 

Solving ( p6| ) for Z = [z, z') with (z, z') = Zj^^) = (z;, z[) we find 

Zo = (0,0), Zp^{zp,Zp), Zoo = (oo,oo), Z± = {z±,z^), (27) 
where Zp is the same as (^, and 
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z± 



l-2g 



^ 1 - (1 - 2g)\e^' ± ^(1 - 25)2(1 - Xe^'^ f + V^— ^-^^ 



(l-AV) 



(28) 



2 (1 - g)\' - (1 - 2g)e^o 

We notice that Zq, Zp, Z^o are also fixed points of the map of the single unit F{z) itself The 
result is already quite interesting. When one of the units is bound to one of the fixed points of the single 
unit the other one is also attracted to the same point in another plane. This can be seen from the formula 



1-3 1 - 2<7 



F{z) 



(29) 
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which follows to (p^. Therefore it is forbidden that, for instance, one falls to the attractor at z = and 
the other falls to z = oo. On the other hand Z± are new and if one of GLM's is attracted to z = z+, the 
other one must falls to z_ , and vice versa. 

In the neighbourhood of a fixed point the behaviour of the mapping is characterized by the nature of 
the Jacobi matrix and its eigenvalues. The general expressions of the Jacobi matrix and the eigenvalues 
associated with our model are given in Appendix. There we also present illustration of the regions of 
parameters which correspond to repulsive, attractive and saddle points. 

We will summarize, in the following, properties of the fixed points. 

Eigenvalues for Zq and Zp are 

2 - Ae** - A'e-*^ 



Ao = Ae'^(l-ff±ff), 



A. 



AA' 



(30) 



They agree with those of uncoupled map when g = 0. 
In the neighbourhood of infinity the eigenvalues are 



A„o = 



(1 - g){a^ + b") ± ^{l~g f{a^~b^Y+Ag^a?b^ 



(31) 



where 



1 „ 1 



lim — . 

2 1 . 1 2;' I — >oc Z 



1 -g + rg 1 - .g ,,.„ 

Notice that the parameter r depends on the direction along which the limits |z| ^ oo and \z'\ — > oo are 
taken. Therefore the behaviour of the map at infinity is not determined uniquely but varies sensitively 
to the direction when there exists the coupling g. When we choose r — ±1, the result is simplified and 
get 



A„ 



\'e-'%l~g±g), if r = 1, 

if r = -1. 



(1-25)2' 

A-t associated with Z± do not have compact expressions but are given by 



(32) 
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A^ 



A 

id 



A'e 



-ie 



l + 2g 



1 - 2g 



1 + 



A'e- 



l + 2.g 



1 - 2.9 



(33) 



is: 



1 



1- AA' 



1 - AA' 



(34) 



2 V(l + A'z+)2 (1 + A'z_)^ 

This result is new since there is no correspondence in the uncoupled map model. 

Given above informations we are ready to draw some pictures about the behaviour of the system under 
a change of the coupling. We observe the following: 

1. Near the fixed points Zq, Zp, Z^o one of two units behaves exactly the same as one of uncoupled 
GLM. The other unit is attracted (or repelled) at the same fixed point. This explaines the coherent 
behaviour often observed in a coupled system. The multiplier of the first unit remains the same 
but the other one changes its scale by the factor 1 ~ 2g. Therefore the second unit increases its 
stability if < 5 < 1 and instability if — 1 < <? < 0. 

2. The fixed points at Z± have no correspondence to the uncoupled GLM. When the coupling is weak, 
(EST) behave as 



z+ = Zp + g 

z- = g 



1 - AA' 



(1 - A'e-»»)2 



1 - A'e 
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(35) 



hence z+ and z_ approach to Zp and respectively. From this we observe that the two units behave 
differently if one falls to z± . 

3. The stability at z± is controlled by |A+| and |A_| respectively. In particular both of them are stable 
as far as |A±| is less than 1. In the weak coupling limit the multipliers A± reduce to 

1 1 



Ah 



2 1 - AA 



; [2 - A'e-"^ - A^A'e^'^ ± (l - Ae"^) (2 - A'e"'^ - AA')] 



which coincide with Ap and Aq evaluated at g — 0. 



B. Linear coupling of ILM 



(36) 



Before studying the coupled GLM, we like to see the effect of the coupling in the case of integrable 
logistic map (ILM). The mapping is given by 



11 



Since the map is linear in z/, it is invertible as long as the matrix cr~^ is regular. Therefore no Julia set 
arises in this coupling. The matrix is singular when g = ^. In this special case the map reduces to 
the single ILM. 

An orbit of the map ( |37| ) can be seen easily if we write as 

a-'^l-2gP, (38) 

2 ^ ^ . , . . ( Z 



and notice that P is a projection matrix satisfying P = P. Every pair of points / I is mapped to the 



z 



plane specified by z + z' = 0. Then the nth map started from I ) yields 



= A"e"'' f ^° ) +A"e™''[(l-25)"-l] ( ) . (39) 

The behaviour of this map is governed by the eigenvalues of the transformation. Corresponding to the 
eigenvalues A± = Ae*®(l — ,g ± g) the eigenvectors of this map are ^j^x) ' ^^'^^ converge either to the 
origin or to infinity, depending whether |A±| < 1 or > 1. These two points are nothing but the fixed 
points of the map of (p7|). 
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IV. JULIA SET OF COUPLED GLM 



We are interested in studying the change of the Juha set when the couphng is introduced. It should 
supply useful information about the effect of the coupling of GLM's. The Julia set in the case of coupled 
GLM, however, must be rather complicated, since two complex variables, hence four real variables, are 
to be considered simultaneously. Moreover a proper definition of the Julia set itself is not clear in the 
case of hyper complex plane. Therefore we must introduce some definition which is appropriate in our 
discussion. For this purpose, in this paper, we propose to call a Julia set an object which is defined 
simply by generalizing z in jl^ to a point in the hyper plane (z, z'). That is, 

the Julia set J{f) of the map ( |2^ ) is a closure of the inverse images 



The set of these points has apparently the following properties, 

1. invariance under the map, 

2. starting from any point of the set the forward map never converges to an attractor. 

The inverse map becomes much more complicated compared with the one of single GLM. In the case 
of the linear coupling in the form of (^6|), however, we can study behaviour of the map analytically to 
some extent. We present in this section how the Julia set of the coupled system will disappear when 
parameters approach to the integrable limit. We also investigate the influence of the coupling to the 
Julia set. Our arguments will be supported by displaying deformation of the Julia set using computer 
simulation. 



Generally speaking it is a difficult problem to find inverse map of a coupled nonlinear system if it 
dedends on more than two variables. In the case of the coupled GLM ([|), however, we can solve this 
problem rather easily and obtain the following compact expression: 




(40) 



of an arbitrary repulsive fixed point 




A. Inverse map of the coupled GLM 



1. General feature of the inverse map of the coupled GLM 




(41) 
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where 

From this expression we see that a pair of points (zi+i, Zj^i) are inversely mapped to four pairs of points. 
We must take into account all of such pairs at every step of the map. 
There are four types of mapping which we represent as 

Here A and B are the operations which map z to A{z) and respectively, according to ([l7|). A single 

operation of the inverse map to a point, say Z — ^ ^/ ^ ? yields the four pairs of points 

M(,)Z, J = 1,2,3,4. (44) 

If we apply this map twice we get 16 points, 

M(,^,,,)Z := M(,^)M(,,)Z, ji,j2 = 1,2,3,4. (45) 

Similarly if we apply the same map n times, we will get 4" points 

^Oi := M^n)Mu,) ■ ■ ■ M^^^^^Z, ji, j2, • • • , j„ = 1, 2, 3,4. (46) 

If we start the map from a repeller Zj, the union of all these points forms the Julia set, which we 
represent as 

J^'''^=lji U M^n..,-..)ZX (47) 

" = \jl,j2,---;jn ) 

Here IJ^^ ^-^ ... ^ means the union of all possible sets of {j\,2i, ■ • • ,in)- 

2. 'Julia set' of the coupled GLM at integrable limit 

It is instructive to study first the case of integrable coupled system, i.e., e = 0. For this purpose we 
write Af(j-)'s in (^^, as 

M(,) =Mf°) J = 1,2,3,4, (48) 

where 



M-:=p., M-:= - _^ a, M- := ^ " U, M- := -Aa, (49) 
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with A being the map which brings an arbitrary point to the fixed point A, and 



^ , E 0\ „ f E 



^(3) - ( 0^ ^)^' %)--(o e)'^- ^^^^ 

In the e ^ hmit, -E(j)'s vanish and the 'Juha set' is simphfied^ 

= U <....:^^\ (51) 

Namely in we replace Ahy p and B by —A. Notice that this result is different from what we expect 
naively as we take the inverse map of ([3^. The latter corresponds to the pa term in (^). The other 
contributions could be found only as we approach to the integrable limit from non-integrable region e 7^ 0. 
This phenomenon was already observed when we studied the single GLM . The situation is much more 
complicated in the present problem because of the existence of the coupling. 

It is not difficult to understand the reason of this happening. In the integrable limit the fixed point 
~ {zp, Zp) moves to (—A, —A) and becomes super repulsive (hence super attractive of the inverse map) 
in the sense that its multiplier turns to infinity. The operator A in M(2) , ^(3) ■, -^^(4) of (^) brings every 
point to this particular pair of points. 

3. Transition of the Julia set to 'Julia set' in the integrable limit 

We now proceed to show that as e approaches to zero a similar phenomena takes place as in the single 
GLM and the Julia set j'^^'^^ converges uniformly to J^3'°\ if the parameters satisfy certain conditions. 
The proof goes as follows. 

Let us assume |Ao| > 1, where Aq is the one of (|3^. This is satisfied if 



, , f 1, if g < 

Then ( ^/ I = ■Z'o is a repeller and on the Julia set. Starting from this point and applying n times the 



inverse map we obtain 4" points of the Julia set M(^j-^ j^ ... j^^Zo to every set of (ji, j2, ■ ■ ■ :in)- 

We are interested in the difference of M(^j^ j^ ... j^)ZQ from M^j^ . ^Zq. The latter is obtained from 
the former by simply erasing Ej^s in (M8) everywhere. The difference, therefore, has to contain E at least 



■'■^This is not a Julia set in the sense of conventional definition Q. It makes sense only by the definition we 
adopted in (|l^. Hence we call this 'Julia set' to distinguish from others. 
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once. We recall that once this function appears in the sequence of the map, all points are mapped inside 
of the disc of radius R^. This also implies that only the last E in the map is sufficient to be considered, 
because the other i?'s in the past are included as informations within the disk which we do not care 
anyway. 

We further notice that A terms in M^^^ 's reset all informations in the past and put to the fixed value 
A. Hence if such term appears in the series of maps, it is exactly cancelled in the difference and we are 
left with only contributions of p terms. Accordingly the difference takes the form 

n 

fc=i 



where we used the notations 



Ci--=P[1 \y ^^--^Pyl 2)' \y C4 = 0. (54) 

Let ( \ and ( \ be one of such pair of points associated with a particular combination of 

(jii ^2, ■ ■ ■ ,jn)- We calculate distances of upper and lower components in the complex planes separately. 
They can be written as 



"^j" ^^j" - 1 ^ ■ ■ ■ 1 '^^Uk ) ^ 

k=l 

k 



(55) 



Each term in the summation over k contributes the factor p" . Let us introduce the notation a^n to 
denote the maximum of the absolute value of the eigenvalues of cr, that is 

:=Max|l, . \ , 1. (56) 



|l-2ff| 

Then we can prove 



< i^^i?. (57) 
1 - \pcrm\ 



if we replace Cj, j = 1, 2, 3 by p 



1 
1 



This result shows that the coefficient of i?^ is finite for all n if \p(J,n \ is less than 1. When 5 < or there 
is no coupling, am = 1 and the condition reduces to the case already studied for the uncoupled map. On 
the other hand, if g is positive, am is larger than 1. In this case the condition becomes |A'| < |1 — 2g\. 
Therefore we can say that the Julia set J^^''^-' converges to J^9,o) uniformly if the parameters satisfy, 
together with (p2|), the following conditions: 

M<{|1-2„|, ito'ls^i- 

A few remarks are in order. 
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1. This condition ( pq ) is not necessary but only sufRcient for the uniform convergence. 

2. In contrast to the uncoupled case the condition |pcr„i| < 1 may not be automatically satisfied when 
|Aoo| < 1- As we see from the expression for Aoo in (|3l| ), it is not easy to signify the condition for 
the infinity to be an attractor of the map. 

In order to see visually how the dependence of the Julia set on the integrability parameter e looks 
like, we will present computer simulation of the inverse mapping for some values of e and fixing other 
parameters. In the case |Ao| > 1, z = is a repulsive fixed point. Hence the Julia set can be generated 
starting from there. Note that in order to visualize it in a finite region of a graph we must perform 
the calculation under the condition \A^\ < 1, so that we do not fail to see any point from our view by 
adjusting the scale. 

Since the Julia set of our system is in two dimensional complex space (z, z'), we are not able to present 
its total figure in one graph. The best we can do is to present its certain projections. Here we will 
show projections of the Julia set to two dimensional real space (i?e(z), Re(z')). The simulations are 
executed at the coupling parameter g = 0.3 and A = 4.0. The following pictures represent the projections 
corresponding to various values of the integrability parameter e. Here the scale of windows is 6 x 6 and 
the center of the windows is at (i?e(z), Re{z')) — (—2, —2). 
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Fig.l Julia set projected on Re{z) — i?e(z')-plane 




e = -0.2 e = 

B. Dependence on the coupling strength 



Now we investigate ttie dependence of tlie Julia set on the coupling strength. As we see in ( pl| ) the 
inverse map depends on the coupling strength g in rather complicated fashion. If we use the variables 
= aZ instead of 2^ = ^^/^, however, the dynamics becomes much simpler. In fact in the {x,y) 
frame the coupling dependence appears in the map only through a given explicitly in (p3[). Wc can 
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employ this frame to analyze the coupling dependence. The effect of the coupling to the (z, z') variables 
can be obtained if we transform the results of the {x, y) variables to those of (2, z') by . 
In order to separate the coupling dependence let us write a as 

cr = l + aP, a , (59) 

1 — 2g 

P is the projection operator defined in (pq). The projection operator P has a peculiar nature. Namely it 



maps any pair of points in the space of two variables ( ^, ) to a symmetric pair of points different only 

z') *° (a(z'-z) 

1. Maps of the form M[j)X 

Starting from the repulsive fixed point Zq, we apply n times the map and see the projection of 
the image on the z plane. There should be 4" points, as given by (p6|), if there is no degeneracy. Before 
going into detail let us first see how the single operations Af(j), j = 1, 2, 3, 4 map an arbitrary coordinate 
X — ^ ^' ^ dimensional complex space ^ 2' ^ ' definition (^ ) we find 

x \ _ ( p{x + ai) + E{x + aO 



by their signs. Hence aP maps I ^, ) to ' '^^^ ^ ^ 



^h^)^ - ( R ) M ) = ( -^-'ku'': ] ' (61) 



A 








A 


A 








B 


B 








A 


B 











A - E{x' - aS) 



where we defined ^ := x — x' . From this result we learn that the operation of B maps all coordinates into 
the neighbourhood of the fixed point —A with radius R^. Since a is only in the function the coupling 
dependence will not be seen unless we watch deep inside of the disk. In this sense the operation of B 
resets the large part of informations in the past. The operation of A, on the other hand, depends linearly 
on a. Multiple operation of A will deform the Julia set depending on the strength of the coupling g. 

2. Maps of the form M(")(Zo) 

In order to see how the Julia set is changed from J*^°'^-', the case of uncoupled map, we study behaviour 
of M("^^Zo in detail. 

If Zk = ( ) is the image of Zq after k times of the map, (|60|) tells us that 



Zk \ _ ( pizk-i + a^fc-i) + Ek-i 



(64) 
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where :— Zj — z'^ and Ej :— E{zj + a^j), E'^ :— E{zj — cn^j)- Solving the recurrence formula 

a = Zk - 4 - p(l + 2a)a-i + Sfc-i - ^Li (65) 

we find 

n-l 

e„ = + 2a)"eo + E /^"(l + 2a)'' {E^^-k-i - , (66) 

A:=0 

from which follow 

n-l 

Zn = p^Zo + ^ p'^En-k-l 
k=0 

-J n n — fc— 1 

+ -p"((l + 2a)"-l)Co + «Ep'- E p'(l + 2a)'(i?„_fe_,_i-£;,U_,_0' 

fe=i ;=o 

n-l 
k=0 

n n — fc— 1 

+ - p" ((1 + 2a)" - 1)^0 - a^p'^- ^ + 2ay {E,,-k-i-i ~ • (67) 

k=i i=a 

The coupling dependence of the first two terms in z„ and enters only through the function E. If we 
ignore the detail inside of the disk of the size i?e , we do not distinguish these terms from the points of 
the Julia set J'^"''^^ of uncoupled system. 

The last two terms both of z„ and z!^ exhibit explicit dependence on the coupling g. The feature of 
the dependence is quite interesting and we will study further detail. The term i p" ((1 + 2a)" — l)^o 
depends on the initial position Zq explicitly. It tends to shift the central position I ^^^^ ) on the 'Julia 



P z'o 

set' j(O'O) = ( jji^i^i j by a small amount as long as the coupling, hence a, is small. The last term is a 
finer correction whose dependence on the initial position comes only through the function E. Using the 
inequality ([l^) the value of this term can be estimated as 

n — k—l 



a 

fc=i ;=o 



2\a\R, /, | 1-H" |l + 2a|"-l \ 

< l-|p(l + 2a)| U"'T^ " |l + 2a|-l j ^''^ 

The value of the right hand side is determined only by the parameters characterizing the map. Therefore 
the information about the initial values of the map is lost in this term. 

It is remarkable that we could get an 'analytic expression' of an orbit after multiple inverse map 
Since A contains a square root function, it is not trivial at all to get any compact expression of an orbit 
for the coupled system. What is meant by 'analytic expression' is that all expressions, but contents of the 
function E, are given explicitly. The maximum value of E is bounded by Re, irrespective to the position 
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of the map. We loose informations inside of the disk at every step of the map. The uncertainty, however, 
does not increase as steps of the map increase, since the difference from the uncoupled case is calculated 
independently at every step of the map. This fact enables us to predict the inverse orbit, as far as detail 
informations inside of the disk are not required. 

3. Maps of general form 

If X is the image of Zq by Af(ji.j2,... j„-s-i) with s between and n — 1, {M^-^-^M(^j'fX; j — 1,2,3,4} 
represent general form of the map Mq-^ j2,...,j„)'^o- We have already discussed the extreme case, s = 
71 — 1, j = 1. From previous experience it is not difficult to see general feature of the map. 

The image of M(^j^X by the map M^^-^ is given by (|6^) if n is replaced by s and Zq is replaced by those 
from @ to (H). 

s-1 

-s,3 



Zn,j = p'Zn-sj + p'^Es-k-l + ((1 + 2a)"' - 1) 
fe=0 

s s — A: — 1 

+ ^Y.p' E p'a + M{Es-k-i-i-E:_,_,_,), 

k=l 1=0 

s-1 

4,, = P^4-s,, + E p'K-k-1 + P' ((1 + - 1) ^„_.,, 

s s—k—1 

-^T.P' E p'{l + 2ay{E,.k-i-i-El_,_i_,), (69) 



where ( , "'''^ ) = M(j}X. We notice the main difference of this expression from the one of (|67|) is that 
the first terms p'^zq and p^z'q in ( |67| ) are points of the Julia set J^^^^ of uncoupled system, whereas 
p'^Zn-s,j and J are not in general. If the last operation of M(^jj was M(4), they must stay in the 

neighbourhood of — /9"^''A, a point in J^^*^. Hence ( ^7'^ ) also remain in the neibourhood of 
Otherwise we do not have full information. 

If there exist points, whose locations relative to those of uncoupled case are not known, the Julia set 
j(9^e) q£ ^Yie coupled system tends to spread compared with J*^"^^^. Detail inspection of the expression 
(|69|), however, will show that such expansion of the Julia set will not continue without limit, as long 
as the parameter p is chosen less than 1 and the coupling strength g is sufficiently small. In fact terms 
including z„_s or Cn~s, hence dependent on the locations, are multiplied by factors of the form p" and/or 
(1 + 20:)" — 1, which tend to push the points to the neighbourhood of the origin. The behaviour of the 
map will change its nature significantly depending if the value of \p{l + 2a) | is greater than 1 or less. 
The critical value is given by 
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The change of the Juha set with the couphng parameter g is represented below. The simulations are 
executed at A = 4 and A' = —0.3. In this case |Aoo| < 1 and |Ao| > 1 hold for small g, hence the Juha set 
exists in a finite region including the origin. This condition, however, is destroyed at some critical value 
of g since the condition of the convergence depends on both |Aoo| and |Ao| and the Julia set may spread 
to infinity. In our choice of the parameters we observe 

1. g= iz^M ~ 0.226 

One of |Aoo(r = — l)|'s becomes 1, hence the Julia set starts to includes the infinity. 

2. g^ gc~ 0.269 

|p(l + 2a) I becomes 1, hence the distance from the uncoupled case fails to be estimated. 

3. 5 = = 0.35 

Another |Aoo(r = — 1)| becomes 1, hence the Julia set tends to expand to infinity. 

4-5=|(l-^)=0.375 

|Ao| becomes 1, hence the origin ceases to be a repeller and the Julia set starts to cover inside of 
the basin of Zp. 

Here we like to visualize the coupling dependence of the Julia set by using the computer simulation of 
the inverse map when other parameters are fixed. Like the previous pictures only projections of the Julia 
set are visible. In the present case we will use the projections to the z plane, so that we can compare the 
influence of the coupling with those of single map. The size of the windows is 8 x 8 and the center of the 
windows is at (i?e(z), Im{z)) = (—2,0) of the complex plane. 



22 



Fig. 2 :Julia set projected on i;-plane 

Fig.2a :g = 0.00 Fig.2b :g = 0.20 




V. DISCUSSIONS 



We have investigated a coupled complex dynamical system which is generated by degree two rational 
map. The subject we considered in this study was the Julia set. It characterizes dynamical feature of 
the system. It is an invariant object in the sense that it does not depend on initial values. It, however, 
changes very sensitively if the parameters of the map are changed. In particular we were interested in 
how the Julia set is influenced when two independent maps are coupled. We have chosen the model of the 
coupling carefully, so that we can study the effect of the coupling, not only numerically, but analytically 
to some extent. The introduction of the coupling tends to spread the Julia set when it is projected to 
the complex plane of the single map. This naturally happens because the projection procedure linearly 
superposes the effects of the other unit of the map. As the coupling strength increases beyond the critical 
value Qc, the Julia set extends to infinity. This was observed also in the numerical calculations. 

Another problem which we were interested in is concerned with the transition of the system from 
chaotic phase to integrable phase. It is a hard task in general to signify such a transition, because 
it requires to analyze some singular functions which appear in the inverse map. When the system is 
complex the inverse functions are seldom to be studied analytically. Therefore it is remarkable that we 
could analyze the transition of the Julia set of our coupled system and could prove the uniform transition 
to the integrable phase. 

Before closing this paper we will discuss briefly an attempt of measuring complexity of the Julia set. 
For this purpose we study fractal dimensions of the Julia set. Since the Julia set is an invariant of the 
map, the fractal dimensions are also invariant once the parameters of the map are fixed. We are interested 
in how the fractal dimensions change as the coupling between two GLM's varies or as the integrability 

parameter e approaches to its critical value. 

There are several kinds of definition called fractal dimension. We will use the box dimension and the 
information dimension in the following, which are relatively easy to calculate from the data of computer 
simulation. 

The box dimension is defined by 

log(jV(^)) 
cLb = hm - — 

5^0 log(l/^) 

where A''(i5) denotes a number of points which constitute the figure under consideration and 5 represents 
the length which decides minimum scale of the measure. 
The information dimension is defined by 

\og{S{5)) 



di = lim 



s^o log(l/5) 
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in which S is the entropy obtained from the probabihty Pj that the mapping across the jth-cell box of 
length S, 

oo 

5(5)^-^P,log(P,). 

In the calculation of the fractal dimensions there arise some problems due to the use of computer. 
Namely, since we can not take the limit 5 — > in the computer, there could remain some ambiguities on 
the values of the calculation. We performed the calculation at different values of S and found that it is 
enough if the value is sufficiently small. To be specific we will fix the minimum scale 5 to be of the 
window length in each computer simulation. In this ratio, we calculate 800,000 points for each plot of 
fractal dimension with a set of parameters. This number of points is determined such that the dimension 
of the random mapping is sufficiently close to 2. 

Another problem which arises in our calculation is due to the fact that the Julia set is an object in 
two dimensional complex, hence four dimensional real space. Very huge number of points are needed 
to calculate the fractal dimensions. In order to reduce the number of calculations we performed the 
calculation only of projections of the Julia set to two dimensional real space. In other words we abandon 
half of the informations about location of each point in the four dimensional space. The values obtained 
in this way do not mean those in the full dimensions. Nevertheless the feature of their variation should 
represent characterization of structure possessed by the Julia sets with different values of parameters. 

The results of our computer calculation of the fractal dimension of the Julia set arc shown in Figs. 3. 
The first two figures a, b are for the fixed values of A = 4 and g = 0.3, whereas the integrability parameter e 
is changed from —1 to 0. They exhibit the gradual decrease of the dimension as the integrability parameter 
e approaches to the integrable limit, hence agree with the analytical result. The same quantities were 
also calculated against the coupling parameter g. 

Figs. 3 c, d present the results of the calculation for various values of the coupling g between and 
0.45, while A and A' are fixed at 4 and —0.3. We observe that both box and information dimensions 
increase as the coupling increases when g is small. Beyond the critical value gc { — in the case of the 
calculation), we observe no significant change of the values of the fractal dimension. This means that the 
expansion of the Julia set to the whole space takes place gradually as the coupling of the system becomes 
strong beyond the critical value. However such an expansion, in turn, causes disappearence of the points 
from our observation. 
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Fig. 3a Box dimension vs the parameter e: 
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Fig. 3b Information dimension vs the parameter e: 
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Fig. 3c Box dimension vs the parameter g: 
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Fig. 3d Information dimension vs the parameter g: 
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Appendix 

The Jacobi matrix of the coupled system is defined by 



d^p{z,z') dip{z,z') 
dz dz' 

dip{z' ,z) d(p{z' ,z) 
dz dz' 



(1 - g)F'{z), gF'iz') 
gF'iz), {l-9)F'iz') 



where 



F [z) = — — = — . . e- 



dz 



(1 + A'z)2 



(71) 



(72) 



The eigenvalues of this matrix are 



1 



A = -(1 - g){F'{z) + F'{z')) ± - v/(l - gfiF'iz) - F'{z')f + Ag^F'{z)F'{z'), 

The behaviour of the dynamical map at infinity is determined by the Jacobi matrix 

d 1 d 1 



9« , 

a 



(73) 



which should be calculated at <^ = i = 0, C = p- = 0. 

In order to get idea about stability of the map near the fixed points, we calculated regions of parameters 
which correspond to repulsive, attractive and saddle points. Under the restrictions of the the parameters 
by A' = A and 9 — 0, we can draw the picture on the (g. A) plane. Figs. a,b,and c show the parameter 
plane associated with the fixed points at Zq, Zp, and Z±, respectively!^. Here the most dark gray region 
represents the parameters corresponding to repulsive fixed points, the gray to attractive and the lightest 
one to saddle points. 
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''''The case of Zaaif = 1) is the same with one of Zq. 
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